Atom focusing by far-detuned and resonant standing wave fields: Thin lens regime. 
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The focusing of two-level atoms in a beam or trap after interacting with both far-detuned and 
resonant standing wave fields in the thin lens and paraxial approximations is considered theoretically. 
The thin lens approximation is discussed quantitatively from a quantum perspective. Exact quantum 
expressions for the Fourier components of the density (that include all spherical aberration) are used 
to study the focusing numerically. The following lens parameters and density profiles are calculated 
as functions of the pulsed field area 6: the position of the focal plane, peak atomic density, atomic 
density pattern at the focus, focal spot size, depth of focus, and background density. The lens 
parameters are compared to asymptotic, analytical results derived from a scalar diffraction theory 
for which spherical aberration is small but non-negligible (9 > 1). Within the diffraction theory 
analytical expressions show that the focused atoms in the far detuned case have an approximately 
constant background density 0.5(1 — 0.635#~ 1//2 ) while the peak density behaves as 3.83# 1//2 , the 
focal distance or time as 6* 1 ( 1 + 1.270 -1 ' 2 ), the focal spot size as 0.744# -3//4 , and the depth of 
focus as 1.910 - 3//2 . Focusing by the resonant standing wave field leads to similar results. However, 
resonant focusing is also accompanied by a new effect, a Rabi-like oscillation of the atom density. 
For the far-detuned lens, chromatic aberration caused by the longitudinal velocity distribution in 
an atom beam is studied quantitatively with the exact Fourier results. Similarly, the degradation 
of the focus that results from angular divergence in beams or thermal velocity distributions in 
traps is studied quantitatively with the exact Fourier method and understood analytically using the 
asymptotic results. Overall, we show that strong thin lens focusing is possible with modest laser 
powers and with atomic beam characteristics that are currently achievable in the laboratory. 
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I. INTRODUCTION 



Evolving from the initial experiments Jljfl], the study 
of atom focusing using standing wave (SW) light fields 
has developed into a broad area of research in atom op- 
tics. A standing wave field acts as a lens as a result of 
a spatially-dependent light shift, focusing atoms into a 
periodic set of lines or dots having widths of the order of 
tens of nanometers and distanced from one another by 
A/2, where A is the wave length of the light. Focusing 
by thick lenses and atomic deposition on a substrate have 
been observed for beams of Na QJ|-01 , metastable He || , 
and Cr The latest achievements in this area are 

summarized in a recent review article JlT| . Perhaps even 
more significantly, lithographic techniques have been de- 
veloped to create nanosurfaces in semiconductors and 
metals using metastable atoms as a pattern template for 
selective etching Q . Similarly, coldjtrapped atoms can 
be subjected to SW light pulses fl^Jbill , creating a peri- 
odic wave packet which will focus along the field propa- 
gation direction at specific times following the atom-field 
interaction. 

Previous experiments with atomic beams have been 
carried out in a thick lens regime for which the atoms 
focus within the laser beam. Since the substrate surface 
is typically placed near the plane of peak intensity of the 



thick SW lens and mechanically fixed to the retroreflect- 
ing mirror, in some sense these experiments are easier 
to set up when compared to thin lens experiments. The 
classical and quantum motion through a thick lens have 
been numerically simulated and compared to experimen- 
tal data in Refs. |^,p"5|-|l8|] . The classical motion and lens 
parameters after thick and thin lenses have also been 
studied in detail within a ray optics formalism p9| . For 
that work the lens characteristics and aberrations beyond 
the parabolic lens approximation were simulated numer- 
ically as were the effects of different atomic velocity dis- 
tributions. However, a comprehensive theoretical study 
of the thin SW lens from a quantum perspective that 
details corrections to the lens parameters as a result of 
spherical aberration, chromatic aberration, and angular 
beam divergence has not appeared previously. 

In this article we stress the flexibility and validity of 
a thin lens approach, where the atoms are detected or 
used for lithography after the interaction with the field. 
We derive exact and computationally expedient results 
for atomic matter waves by quantizing the center-of-mass 
motion from the outset, a process that naturally accounts 
for spherical aberration. Our hope is to stimulate inter- 
est in thin lens experiments by introducing a straightfor- 
ward, yet rigorous theoretical framework for the problem 
and by deriving analytical results for lens parameters to 
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facilitate lithographic and atom optical configurations. 
The thin SW lens creates high contrast periodic struc- 
tures in free flight and opens up numerous possibilities 
for Fourier atom optics to image and manipulate mod- 
ulated matter waves. This approach is highly desirable, 
especially when coupled with advances in Bose-Einstein 
condensate (BEC) atom laser sources [ fl4|]2"o| with large 
fluxes and narrow velocity distributions. To our knowl- 
edge similar studies have been carried out only numeri- 
cally using a Fresnel diffraction theory in the coordinate 
representation for a "doughnut" -mode optical field pj| , 
standing wave field |2^] , and conical lens |2^] . 

Rather than using the diffraction theory simply for nu- 
merical propagation of the matter waves, we extend its 
usefulness by deriving approximate analytical expressions 
for the SW field lens parameters. The Fresnel-Kirchhoff 
integral, incorporating the lowest-order spherical aberra- 
tion, gives a relatively simple form for the atomic wave 
function near the focal plane and SW intensity extrema. 
Asymptotically for large atom-field pulse areas, this inte- 
gral leads to a universal atomic density profile, revealing 
the behavior of the focal distance, peak atomic density, 
focal spot size [half-width at half-maximum (HWHM) 
of the atom density at the focus], depth of focus, and 
background density as a function of the pulse area. The 
analytical asymptotic results can then be compared to ei- 
ther numerical asymptotic approximations of the Fresnel 
integral or to the exact atomic distributions and lens pa- 
rameters, most easily evaluated using a Fourier approach. 

The time evolution of a free particle's wave function, 
when subject to a spatially periodic initial condition such 
as that imposed on an atom beam after passing through 
a SW field or microfabricated grating |l3|j2^ , |25| , can be 
expressed in terms of a Fourier expansion in the (plane 
wave) spatial harmonics of the modulation period. As 
such, the results we derive are closely related to the canon 
of research on periodic light optics |2(| . In the thin lens 
approximation the Fourier components can be calculated 
exactly for the wave function and density and then super- 
posed to evaluate the atomic spatial distribution: for this 
approach the calculations effectively take place in a mo- 
mentum representation, whereas a numerical integration 
over the free particle propagator gives results in a coor- 
dinate representation. The expressions for the density's 
Fourier components that we use below were obtained re- 
cently for far-detuned j27| and resonant |28j SW fields. 
To our knowledge this exact solution has not been used 
previously to examine periodic focusing parameters, nor 
to compare with approximate solutions of focused waves. 
However, the Fourier components and the degradation of 
this type of focusing have been proposed recently as a 
way to study the onset of many-body effects in conden- 
sates pg| |. 

The Fourier approach has several advantages. First, 
the time dependence of the Fourier coefficients directly 
reveals the quantum nature of the periodic atomic den- 
sity as it varies explicitly with the atomic recoil energy 
of integer pairs of photons. Second, summation of the 



Fourier series leads to fast convergence of the Fresnel so- 
lutions to arbitrary accuracy for all times when compared 
to the numerical integration of the full Fresnel-Kirchhoff 
integral. With this Fourier method for example, we can 
reproduce the main Fresnel diffraction results of Ref. |22] ] 
without having to integrate the partial differential equa- 
tion or equivalent integral equation on a lattice. Fur- 
thermore, the Fourier method can be employed for the 
general case of typical atomic beam experiments, where 
the incident beam or trapped atoms are most accurately 
specified by a density matrix in momentum space, such 
as a thermal velocity distribution. The implications are 
that the spatial extent of the atomic distribution is much 
larger than the wavelength of the standing wave and the 
wave function is not transform limited (except perhaps 
for the case of a condensate), so that the focal densities 
are conveniently found as a function of this initial density 
matrix. As a direct result, while the failure to account 
for finite beam size effects with the Fourier series limits 
our calculations to the near-field and Fresnel diffraction 
regimes, it is relatively unimportant for atom focusing 
immediately following the SW pulse. We can then calcu- 
late the effects of chromatic aberration and angular beam 
divergence by averaging the exact Fourier solution over 
any longitudinal and transverse atomic velocity distribu- 
tions. 

In addition to lithographic detection schemes, one can 
detect atomic distributions modulated by SW fields by 
backscattering a pulsed traveling wave probe off of the 
density either (a) at a certain distance from the grat- 
ing for the beam case or (b) after a certain time for the 
atoms released from a trap. The scattered signal is sen- 
sitive to the lowest-order Fourier component of the 
modulated density. Previously, the backscattered signal 
from a cold vapor subjected to strong SW pulses was ob- 
served in Rb in a ground state echo experiment Q . The 
excellent agreement between the Fourier theory and ex- 
periment that was achieved in that case insures that our 
theoretical approach can be applied to the atom focusing 
problem as well. 

A new regime of atom focusing arises if the SW field 
is resonant with the transition between a ground and ex- 
cited state and the time of interaction is less than the 
excited state lifetime. After propagating from the SW 
field for a time just long enough for the excited state 
to decav, the atom spatial distribution consists of two 
parts the stimulated density, caused by direct am- 
plitude modulation of the ground state wave function by 
the field, and the spontaneous density, produced by exci- 
tation to and decay from the excited state. As a result of 
the periodic Talbot revival, or self-imaging 26 3l] , |33| , |3^] , 
of the atomic spatial distribution, free evolution of the 
stimulated part leads to atom focusing not only at the 
focal distance Lf, but also at the distances j'Lt/2 ± Lf, 
where j is a positive integer, Lt — A 2 /2XdB is a Talbot 
length, and X c ib is the atom de Broglie wavelength. In 
contrast to the stimulated part, the spontaneous modula- 
tion decays at a distance on the order of Lt as a result of 
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the Doppler shift associated with spontaneous emission. 
At a sufficiently large distance, the modulated sponta- 
neous part disappears, and one can observe and analyze 
atom focusing by a resonant SW field. This analysis then 
reveals a new effect, a Rabi-like oscillation of the focal 
parameters with the atom-field interaction strength that 
arises from an interference between different components 
of the matter wave. 

This article is arranged as follows. In Section II we 
consider the focusing by a far-detuned field acting as a 
phase grating. Focusing by a resonant field acting as an 
amplitude grating is studied in Section III. Section IV is 
devoted to the influence of chromatic aberration (longitu- 
dinal velocity distributions) and angular beam divergence 
(transverse velocity distributions). The quantitative re- 
sults and illustrative examples are discussed in Section 
V. 



II. FAR-DETUNED STANDING WAVE LENS 

Atom optical experiments can operate in the spatial 
or time domain. For example, in the time domain a 
vapor of cold atoms interacts with one or more radia- 
tion pulses forming spatial gratings in the x direction. 
Correspondingly, in the spatial domain an atom beam 
traverses one or more optical elements or interaction re- 
gions. For a monovelocity beam propagating along the 
z axis with longitudinal velocity U and with the optical 
elements aligned in the x direction, the spatial domain 
configurations can be analyzed in the time domain if cal- 
culations are performed in the atomic rest frame moving 
with velocity U . In this frame the optical elements ap- 
pear as interaction pulses. As a result, our calculations 
are restricted to the time domain without loss of general- 
ity (and are adapted to account for a longitudinal veloc- 
ity distribution in Section IV). We consider the focusing 
effects of a single pulse. 

The atom optical elements couple to the center-of-mass 
degrees of freedom of an atom with mass M and de 
Broglie wavelength \ dB = 2nh/(MU). When the SW 
laser field, 



is the two-state wave function for the ground (i/j(x)) 
and excited (ip e (x)) states, p x is the center-of-mass mo- 
mentum operator, lu is the atomic transition frequency, 
f2, k = q/2 = 2ir/\, and g (t) are the frequency, propaga- 
tion wave vector, and pulse envelope function (centered 
at t = 0, having peak value of unity and duration on 
the order of r), respectively, of the field, (7+(er_) is the 
atomic raising (lowering) operator associated with the 
transition \g) — > |e) (|e) — > \g)), and x = —^E/2h is a 
Rabi frequency for the \g) — ► |e) transition with dipolc 
matrix element [i. We have assumed for simplicity that 
the Rabi frequency x is real. 

In the spatial domain this Hamiltonian is written in the 
atomic rest frame for the slowly- varying wave function 
in the paraxial optics limit, 



MU » (p x ) , (pi) - (p x f 



(4) 



for t > the distance from the interaction region in the 
lab frame is L = Ut. In the time domain a paraxial ap- 
proximation is unnecessary, and is simply the wave 
function. Experimentally ||, for focusing it has proven 
advantageous to optically pump the atoms into an initial 
ground state with magnetic quantum number rrij = ±j 
and use circularly polarized fields to avoid multiple Rabi 
frequencies in the state dynamics, justifying our two-level 
approximation. 

Pure phase modulation of the atomic ground state dur- 
ing the interaction occurs in the far detuned case 



A > max (r, t~ 



(5) 



where A = fl — u is the atom-field detuning, and T is the 
excited state decay rate. Often, a steady-state light shift 
potential, derived from the Hamiltonian (g), has been 
used as an effective atom- field interaction [135(1 < 
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8| X |V(*W (kx) 
(r/2) 2 + A 2 



(6) 



E (x, t) = Ee~ int g(t) cos (kx) + c.c. 



(1) 



While this potential is important for smaller detunings 
and/or larger intensities, it is strictly valid only for 
rr > 1. On the other hand, the dressed state poten- 
tial, 



drives the atomic transition between the ground state 
\g) and the excited state |e), the Schrodinger equation, 
ihd^/dt = H^, governs the motion. The Hamiltonian 
in the rotating-wave approximation is 



H = 



9 

Px 

2M 



TlLJ 



where 



+2h X g(t) cos (kx) (e~ mt a + + e int a 



( i/) e (x,t) 
U(x,t) 



(2) 



(3) 



^V 1 - U > \\ S >ril)^(l>:r)/^. 



(7) 



neglects spontaneous emission and assumes that the 
atom in its ground state adiabatically evolves into one 
of the two dressed states. For our work we follow the ex- 
perimental findings of Natarajan and coworkers || . They 
show that focusing is improved by taking a short inter- 
action time and a large detuning relative to the Rabi fre- 
quency and decay rate. In our case these conditions avoid 
spontaneous emission (diffusive aberration ]l7| ) both 
during and after the atom-field interaction. In this limit 
the two potentials, Eqs. (g) and (0), reduce to the same 
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effective potential, ignoring spatially independent energy 
terms. The Hamiltonian in a field interaction represen- 
tation for the ground state wave function ip(x,t) after 
adiabatically eliminating the excited state for A 3> 2 |%| 
and 4IY|x| 2 /A 2 < 1 (ip e (x,t) ~ 2 X /A ~ 0) is 



H = 



2M A 



-g 2 (t) cos(gx). 



(8) 



If the incident atom wave function is uniform 
[ip (x,t = 0~) = 1] , then just after a single interaction 
the wave function is given by 



tp(x,t = 0+) = exp [i (6»/2) cos (qx)] 

OO 



where 



4jxl 
A 



dtg 2 (t) 



(9a) 
(9b) 

(10) 



is an effective pulse area for the far-detuned atom-field 
interaction, and J n is a Bessel function of order n. Equa- 
tion ( pa| ) is valid in the Raman-Nath |36| or thin lens 
approximation, for which the p 2 /(2M) term is ignored 
during the interaction so that the field acts as a standing 
wave phase grating for the atoms. A standard condition 
given for the Raman-Nath approximation for thin lens 
focusing is p7p2| 



(11) 



\e\u g r/2 < i. 

The two-photon recoil frequency is 



kg 2 
2M 



2tt- 



U ' 



(12) 



and Lt = A 2 /(2Ads) is the Talbot distance for the atom 
beam. The momentum eigenstates, exp[ingx], which are 
coherently superposed to form this wave function, each 
have the free particle energy n 2 Tiuj q . 



A. Thin versus thick lens regimes 

Before we proceed, a clarification is needed to empha- 
size the differences between the thin lens (Raman-Nath) 
and thick lens regimes of the SW atom lens. Some of the 
quantitative differences have been examined by Hcnkcl 
and coworkers Kfl] in relation to Fraunhofer diffraction 
of atoms and from a ray optics (classical) perspective by 
McClelland for focusing |L9|] . For focusing considerations 
we assume square pulses, g(t) = 1 for — t/2 < t < r/2 
and zero otherwise, in order to obtain quantitative re- 
sults. We also require 



Conditions ( |l3| ) are necessary for high contrast, thin lens 
atom focusing, where we are interested in the atoms af- 
ter propagating through the lens. (The pulse shape will 
affect the coefficients of our results, not the scaling with 
9 and w q T.) 



The Raman-Nath regime leading to condition (11) is 
normally defined as an interaction for which the aver- 
age kinetic energy gained by the atoms remains much 
smaller than the interaction strength coupling the mo- 
mentum components, (p 2 '(t) / '(2M)\ <C ?i|x| 2 /|A| in the 
square pulse case. If this condition is violated while 
( |l3| ) holds, the lens is thick, and the atoms can fo- 
cus within the interaction region near the time t ~ 
-t/2 + Tr(\x\ 2 ui q /\A\)-^ 2 /4 < t/2. If the atoms have 
not focused completely by the end of the pulse, they will 
exit the interaction region amplitude modulated. To low- 
est order in 9ui q T, the amplitude correction which multi- 
plies the wave function ( |9a|) immediately after the square 
pulse can be calculated to be exp[9u q T cos (qx) /4]. 

However, the atoms also acquire an additional, 
spatially-modulated phase shift, 9 2 u q T cos (2qx) /24 to 
lowest order in 9 2 u q T. This is called the WKB correction 
by Henkel and coworkers |5?J and has its semiclassical 
origin in the harmonic motion of the atoms during the 
interaction. From a quantum perspective these ampli- 
tude and phase changes are caused simply by the kinetic 
energy acquired during the interaction. Heuristically, us- 
ing the wave function which is evolving in the SW field, 
if) (x, t) = exp[— i2 |x| 2 t cos (qx) /A], the correction which 
multiplies ip (x, t/2) can be written as 



-i l r/2 p 2 
exp[— / dtip* (x,t) tt-tjV (x,t)] 



-r/2 



2M 



(14) 



exp [6iJ q T COS (qx) /4~iuj q T0 2 (l - cos(2gx))/24], (15) 

giving the rigorously correct result. If we are interested 
in far-field diffraction, the momentum state wave func- 
tion can change significantly from the spectrum of Eq. 
@ if 9 2 Lu q T/24 > 1 even if condition (|ll|) is satisfied. 

In the Fresnel focal region of interest here, while the 
spatially-dependent phase shift is crucial in determin- 
ing the thin lens properties of the SW field, we require 
\0\ /2 ^> 9 2 uj q T/6 to assure the dominance of the Raman- 
Nath wave function (|9a|) near the focus. In other words 
a corrected Raman-Nath condition, 



|0|ovr/3« 1, 



(16) 



LU q T < 1, |0| > 1. 



is sufficient to observe thin lens focusing, but the addi- 
tional thin lens condition, 2 ui q T cos (2qx) /2A<^ 1, may 
be necessary to ignore corrections to Eq. (|9a|) for other 
observables, like the far-field diffraction pattern or the 
time-dependent behavior of Fourier components of the 
density in echo configurations. Putting the results to- 
gether, the corrected wave function to lowest order in 



(13) 9 2 uj q T, 9u) q T, 
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t/j (x, r/2) = cxp [i (6/2) cos (qx)} (17) 
x exp [i8 2 uj q T cos (2qx) /24] exp [9oj q T cos (qx) /4] , (18) 

could be incorporated into the theoretical work below 
if necessary. We have verified these quantitative re- 
sults using a Crank-Nicholson technique to integrate the 
Schrodinger equation for the far-detuned SW Hamilto- 
nian (§f) numerically on a lattice |3q] . 



B. Results for far-detuned, thin lens focusing 

Returning to Eqs. (^) as the initial condition for the 
free motion, since only the wave function phase has been 
changed during the interaction, the total atom density, 



P (x,t) = \lp{x,t)f + |Ve (M)| - \ll>{x,t)Y 



(19) 



is initially uniform. As different Fourier components of 
the wave function acquire different phase shifts, <p n {t) = 
n 2 u) q t, during the free evolution, the atom density for 
t > becomes spatially modulated. The period of the 
spatial modulation is equal to A/2 = 2n/q. Transferring 
to the dimensionless variables, 



x — ► qx, t — > U) q t, 

for t > one finds l2(i|J27j| 



ip(x,t) = i n Jn[0/2]e l ^ nx - n2t \ 

n— — OO 
OO 

p(x,t) = J ™ [Osm(nt)]e mx . 



(20) 

(21a) 
(21b) 



n— — oo 



This result, exact in the thin lens approximation, is used 
below for the numerical study of atom focusing. We re- 
fer to the exact calculations by this Fourier technique as 
Method 1 in the text and figures that follow. 

Arbitrarily precise values of the focal time tf (giving 
the focal plane position for the beam, Lf = Utf) for a 
given 9 are defined by the first maximum of the d ensit y 
along x — 0, p (0, t), as a function of 9 using Eq. ( 21b ). 
Given the focal time, we can further characterize the lens 
by the density profile at the focus p(x,tf) and its peak 
p(0,tf), the spot size of the focus w, the depth of fo- 
cus At , and the background density p(w,tf) |59|. The 
spot size w (HWHM of the density profile p(x,tf), where 
p(0,tf) is the maximum) can be defined implicitly as the 
smallest positive root of the equation 



p(w,t f )= x / 2 p(0,t f ) 



(22) 



Similarly, the dept h of f ocus (confocal parameter) can be 
defined using Eq. (|21b|) as the time window, 



(23) 



within which the density along x = rises from half its 
peak value at f_ to its peak value at tf and back again 
at t+, given implicitly by 



P%t±)= 1 /2p(0,tf). 



(24) 



This region is not symmetric with respect to the focal 
time tf as a result of the spherical aberration of the lens. 
(Note that the choice of 9 as positive is unimportant even 
though it implies a red detuning of the field. The 9 < 
case is identical but shifted in x by it.) 

Several of these exact lens parameters, as calculated 
by Method 1, are compared to the approximate diffrac- 
tion theory below. Equations (|21b|) will be modified to 
account for finite beam divergence and chromatic aberra- 
tion in Sec. IV; by its nature the Fourier method includes 
spherical aberration (anharmonicity in the SW potential) 
to all orders in the lens curvature. Although Eqs. j2l| ) 
are exact, they offer no transparent possibilities for ob- 
taining analytical forms for the lens parameters as func- 
tions of 9. In contrast, a diffraction theory can be used 
to find approximate, asymptotic expressions (9 ^> 1) for 
the ^-dependences of the lens parameters. 

Rewriting Eq. ( |21a| ) as an integral in the coordinate 
representation, we can express the wave function as 



ip(x,t) 



dx'G(x-x',t-t')ip(x',t'), (25) 



where the Fresnel-Kirchhoff propagator of the free atom 
motion is given in the dimensionless variables (^p| ) by 



G(x,t) 



(AttU) ^ exp (ix 2 



/At). 



(26) 



The integrand wave function is taken as ip (x',t' = 0) = 
exp[i(9/2) cos(x')]. In this form we are first interested 
in the wave function behavior near the focal points, 
x m = 2irni, for integer m and times t immediately follow- 
ing the interaction. The relevant time scale for focusing 
will become apparent shortly. When 9^-1, the main 
contribution to the integral ( pa ) at these times 

arises from small values of \x 1 ~ — x m \, where the poten- 
tial in Eq. (||) is nearly harmonic. We choose to expand 
around the point x' = (m = 0) and downplay the pe- 
riodicity of the wave function for the diffraction theory. 



Replacing cos (x') by 1 
center, 



"/2, the atom density at the 



p(o,t)~(i-tey 



(27) 



contains a singularity at t — 9~ 1 which determines the 
approximate position of the focal plane. 

To obtain a finite value for the atomic density, spher- 
ical aberration (i.e., anharmonic terms in the potential) 
must be considered. Expanding cos (x') to the x term, 



At = t + - t- for t- <t f <t 



omitting the phase factor exp (t 

scaled position £ = 1 /2(9/3) ^ 
able, we find that the asymptotic wave function in the 
vicinity of x = is 



2 /At) , and choosing the 
x 1 as an integration vari- 
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ip (x, t) ~ (3/6») 1/4 {-Kit)' 1/2 exp (iG/2) f [x, Q) , (28a) 

/CO 
d£ exp (-ix£ + iu>£ 2 + i^) , (28b) 
-OO 

x = {3/9) lfi x/t, (28c) 
Cv={3/9) 1/2 (t- 1 -9). (28d) 

The main contribution to the wave function comes 
from the region where the integrand phase originat- 
ing from the anharmonic term is of the order of unity 

UCmaxI ~ 1 ° r l^maxl ~ ®~ • To determine the valid- 
ity of Eqs. (pq), the next anharmonic term (~ x m ) must 
produce a small addition to the integrand phase in this 
region, implying fe^ ax < 1, or 



I- 1 /* 



< 1. 



(29) 



This condition defines the proper asymptotic limit - a 
wave function correction of relative weight 6~ / 2 in Eq. 
(28a) could be included to increase the calculation's ac- 
curacy. We assume condition ( p9] ) holds. 

>From Eqs. ( p8|) the time-dependent density, 



p(x,t) 
F{x,u>) 



\^(x,t)\ 2 r 
l/(^)| 2 



(3/9) /2 (irty 1 F(x,Q) (30a) 

(30b) 



contains both the slowly- varying t^ 1 dependence and the 
sharp dependence of F (x,uj) in the vicinity of t = 9~ l . 
At the focal center (x = 0), the integral ( [28b] ) can be 
expressed analytically (see Ref. No. 3.696) through 
fractional Bessel functions as 



/(o,a) 



3 , 1 1 

/2 7r \u)\ 12 exp 

3 / 2 



7T)/8] 



J_ 1/4 (^ 2 /8) + l '/ 2 si 5 n (w) Ji /4 ((i 2 /8) 



(31) 



giving an analytical asymptotic expression for the density 
as a function of time along x = 0, 



P (o,i)~(3/0) 1/2 (Trt)- 1 !/^*)! 5 



(32) 



The lens parameters in this asymptotic diffraction the- 
ory are defined in the same way as they were for the exact 
Fourier theory above. They are found graphically or nu- 
merically by evaluating Eqs. @, @), and @. The 
focal time tt for a given 9 (which is not given by t = 9~ v ) 
can be found as the first maximum of p(0,t) from Eq. 
(j32|). Using this focal time tf, we further characterize 
the SW lens by the peak density p(0,tf) evaluating Eq. 
([32"|), by the density profile at the focus p(x 7 tf) found 
from Eqs. (§|) and ©, by the spot size w (HWHM of 
the density profile p(x , tf)), and by the depth of focus At 
found from Eqs. (30a) and (^). 

Calculations performed with these asymptotic diffrac- 
tion results, only restricting x to be — 27rm, are 
referred to as Method 2. In addition to the Method 2 
asymptotic results, we can derive expressions which are 



valid both for x ~ x m and for the time restricted to 
be near the focus, t ~ 9~ 1 , namely within the depth of 
focus. This further approximation, which is referred to 
as Method 3, leads to analytical expressions for the lens 
parameters in the asymptotic limit. 

In order to proceed, we need to determine the peak of 
the function F (0,a>) jt from Eq. ( |30a| ) near t ~ 9~ 1 for 
9 / 2 ^> 1. First, times within the depth of focus \t — 
are roughly determined by the requirement \ui\ ~ 1 or 



\t 



(33) 



and scale as ~ 9 ^ 2 times the focal plane position. For 
6* 1 / 2 > 1, Eq. @ justifies the inequality 



< t. 



(34) 



Therefore, it is sufficient to replac e th e slow ly-v arying 
time dependence i" 1 by 9 in Eqs . (28c) and (30a), tak- 
ing the lowest order limit of ( [28d| ) near t ~ Q~ x . In other 
words we replace Eqs. (28c) and fl28cj ) by 



x a ~ 3 fi 9 3/i x and (35a) 
Co a - 3 1/2 6> 3 / 2 (9- 1 -t) , (35b) 

respectively, to arrive at the density near the focus, 

p(i,t)~3 1 /» 7 r-V/ a F[x ,w ]. (36) 

The subscript a reminds us these results are asymptotic. 
Equation ( p6[) can be evaluated numerically. Unlike Eqs. 
( |30a| ) and (|32|), Eq. ( |36| ) is independent of an explicit 
dependence on t and, therefore, is a uni vers al (scaled) 
density function. Inserting Eqs. (|3ll ) and ( 35b ) into (|3€|), 
one finds that the asymptotic time evolution of the den- 
sity along x = can be written analytically as 



p(0,t)~ 
x J_i /4 



)-3 



(s^) 1 / 2 ^!^! 



(37) 
(38) 



and is shown in Fig. [l]. 

Approximate expressions for the focal parameters 
are derived by finding the maximum of the function 
F (0, uj a ) = 1/(0, u) a ) | 2 either numerically using Eq. ( |37| ) 
or graphically from Fig. [ij We find that the maxi- 
mum F(0,ujf) « 6.94 occurs at a)/ » —2.20. He nce, the 
asymptotic focal time is given by inverting Eq. (|35b|), 



1 - Co f (30) 



-V2 



1 + 1.270" 



-V2 



(39) 



The asymptotic atom distribution at the focal plane, 
p(x,tf), putting Co a = Cof into Eq. (|36|), is also plot- 
ted in Fig. |l|. For the peak density we have 



p{Q,t f ) - 3 1/2 TT- 1 F(0,Co f )9 1/2 ~ 3.836» 1/2 . 



(40) 



The asymptotic focal spot size w can be expressed 
through the half-width if of the function F (x a ,ujf) , 



6 



which is found to be if 
gives the expression 



0.979. Inverting Eq. (|35a 



w 



3 A 



0.7446l~ 



3 A 



(41) 



The asymptotic depth of focus, At, defined by Eqs. ( |23| ) 
and (p4]), is evaluated using Eq. (p7|) by finding the 
two values 5)± on either side of tD/ « —2.20 for which 
F(0,Q ± ) = F(0,u) f ) /2 w 3.47. These are Cj+ w -3.42 
and u)_ w —0.115, giving a depth of focus from inverting 



Eq. ( |35b| ) of 

At 3- 1/2 (w_ 



r 3 /2 ~ i.9i6»- 3 /2_ (42) 



Equations (J35|)-p2|) constitute the Method 3 results de- 
scribing the SW lens near x — 2ttiti and t = 0~ 1 . 



© 

a 

■ — - 

CD 



o 

II 

es 

>- 
H 

HH 

GO 

z 
w 

Q 









- « 




\ (a) 




- 








- N 


t r 






- c 

- 


LO - 5 


5 10 




TIME, 3 1/2 9 3/2 


(t - e-0 



Q_ 

CD - ve 

H 

GO 

z 

— 

•< 
u 
c 







(b) 


2w M 








FIG. 1. Scaled asymptotic (as 8 — » oo) atom densities in 
far-detuned focusing, (a) Time evolution of the atom den- 
sity at x = 0, p(0,t). (b) Density profile at the focal plane, 
p(x,tf). Time zero in the graph corresponds to the focal point 
in the harmonic approximation. 

In much the same way, we can use Methods 2 and 3 
to derive expressions for the background density, defined 
as the density p{xb,t) at the midpoints between periodic 
focuses, Xb = x m + ir. We need to find the wave function, 
Eq. (p5j), at the point Xb — it. Since no singularity arises 
in the integral for t < tf when the integrand is expanded 



around x' = n, only the harmonic term is needed to find 
the wave function, where cos (x') ~ (x' — 7r) 2 /2 — 1. We 
find the approximate wave function, ~ (1 + t9)~ 1 /' 2 , and 
background density, 



P (n,t) ~ (1 + tey 1 , 



(43) 



correct to order 8~ 2 . This is the Method 2 result. This 
expression allows us to find the asymptotic density con- 
trast of importance for lithography, c(t). The contrast 
is defined to be the ratio of the atomic density at the 
focal points x m to the background, 



c(t)=p(0,t) /p(n,t) 



(44) 



where p (0, t) is given by Eqs. @ and @ for Method 
2. 

In addition, putting the asymptotic focal time (|39| ) 
from Method 3 into Eq. (fl3]), the asymptotic background 
density at the focus is p{-K,tf) ~ (1 + tfd)~ 1 , or 



p(n,t f ) - 0.5(1 - 0.63561" 1/a ). 



(45) 



^From Eq. ([hi]) the asymptotic contrast ratio at the focus 
is (within the accuracy of this calculation) 



1 pfatj) 



Equations (flit) and (Eq) are a Method 3 result. 



(46) 
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FIG. 2. Far-detuned focusing. Comparison between ex- 
act densities by the Fourier method (Method 1, solid lines) 
and approximate asymptotic densities near x — (Method 2, 
dashed lines) for different values of the pulse area 9. (a) Time 
evolution of the atom density at x = 0, p(0,t). (b) Density 
profile at the focal plane, p(x,tf). 



7 



To summarize, we now have three methods to find the 
lens parameters, each with different levels of numerical 
effort required to use it. The relevant equations for these 
methods are listed in Table 1: 

Table 1. Techniques and equations for three different 
methods of calculating thin SW lens parameters. The 
numbers in the Table correspond to equations in the text. 



Method 


Numerical 
Sum 


Numerical 
Integral 


Analytical 
Expression 


1. Exact Fourier 


M 






Asymptotic 
2. \x - mir\ < 1 
t > 




1 


| 


1 
1 


Asymptotic 
3. \x — mir\ <C 1 
t-6- 1 




E 


|5 


M 

4| 




PULSE AREA 9 



fa 

u 
< 

H 



— 

< 

o 

fa 



© 
d 
d 
d 



1(b) 




A V 

A 

w 
w 




_ ou 





10 



20 



30 



40 



o 

w 

H 

O 

*J 
O 

n 



FIG. 3. Pulse area dependences of the lens parameters in 
far-detuned focusing: (a) Peak focal densities p(0,tf) (curves 
p) and spot sizes w (HWHM of the density profiles at the focal 
planes). (6) Focal distances tf and depths of focus At. Com- 
putational methods: 1. Exact Fourier (solid), 2. Asymptotic 
diffraction near x — (dashed), 3. Asymptotic diffraction 
near x = and t = d' 1 (dot-dashed). Note that the three 
methods give almost identical result for tf. 

In Fig. U the exact atom density at the center p (0, t) 
and atom distribution at the focal plane p(x,tf), as cal- 
culated by Method 1, are compared with the approxi- 
mate expressions from Method 2. The convergence of the 
Method 2 result as 9 increases is evident. In Figs. || the 
pulse area dependences of the focal parameters, as calcu- 
lated by the three methods, are compared. The accuracy 
of Method 2 is generally better than that of Method 3, 
but all methods converge for 6* 1 / 2 3> 1. These results 
are discussed further in Sec. V, including the easily ex- 
plained, constant offsets of the asymptotic peak densities 
in Fig. |a. 



III. RESONANT FOCUSING 

We now consider the density modulation and focusing 
caused by a resonant standing wave field acting on a two- 
level atom. If the pulse duration satisfies the inequality 



(47) 



[the field area 9 for a resonant SW field is defined below 
by Eq. (|5C|)] , then during the interaction the ground and 
excited state wave functions, ip (x, t) and ip e (x, i), evolve 
in an interaction representation according to 

ip(x,t) = — 2i\g (t) cos (kx) tp e (x, t) (48a) 
ip e (x,t) = -2ixg(t) cos(fcc) ij) (x,t) . (48b) 

By satisfying Eq. (|47|) , we are assured that the reso- 
nant atom-field interaction is in the Raman-Nath regime 
and that the pulse duration is shorter than the excited 
state lifetime, avoiding the complications of saturation 
and momentum space diffusion. If before the interaction 
the incident wave function is uniform in the ground state, 
ip (x,0~) = 1, then just after the interaction 



V>eO,0+) 



-zsinp/2) cos(kx)} 
cos[(0/2) cos (kx)] 



(49) 



where the effective pulse area for the resonant atom-field 
interaction is defined as 



(50) 



9 = A X dtg (t) 



PULSE AREA Q 



The subsequent free space evolution and radiative de- 
cay of the system were recently analyzed |28| using the 
Fourier method for the closed two-level scheme (i.e., the 
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excited state |e) decays only to \g)) O]. However, the 
results that are relevant to the problem discussed here 
pertain to any excited state decay scheme. In the follow- 
ing paragraphs we summarize these results in order to 
apply the three methods of Sec. II to the resonant case. 
Our goal is to separate two terms in the free evolution of 
the ground state: (1) the strict Hamiltonian evolution, 
resulting from the ground state amplitude produced co- 
herently by the SW pulse, and (2) the terms which re- 
sult from spontaneous decay of the excited state. Such 
a separation is straightforward by the Fourier technique, 
Method 1. Once this is done, we can use a diffraction 
analysis to obtain asymptotic results for the focal den- 
sity caused by the Hamiltonian term, ip(x, + ) in Eq. 
(]49|). The diffraction analysis, Methods 2 and 3 from 
above, leads to analytic expressions for the resonant SW 
lens parameters and, anticipating the results, sheds light 
on the Rabi-like oscillation of the lens parameters as a 
function of 8. 

After the excited state decays (in a time t > T _1 fol- 
lowing the pulse), the total atomic density p(x,t) in- 
volves only the ground state and is a periodic function of 
x having period A/2. The atomic system can no longer be 
described simply by a wave function since we must trace 
the full density matrix over the spontaneously emitted 
photons to form the reduced density matrix p(x,x',t). 
However, the part of the density matrix which was pro- 
duced coherently by the pulse can still be described by 
evolving the wave function ij)(x, + ) in Eq. ( ]49| ) into 
tp(x,t). This was done in Ref. P8| . The total den- 
sity p(x,t), which is formally the diagonal component 
p(x,x,t) of the reduced density matrix, can be expanded 
in a set of Fourier components, 



(D) 



n— — oo 

= p^ (x,t)+pW (x,t) 

OO 

p^(x,t) = £ pi s Ht)e in ** 



inqx 



n— — oo 
oo 



fat) = E 



(51a) 
(51b) 

(51c) 



n— — oo 



Each ground state Fourier component, p n (t) — p n (t) + 
pn (t) i consists of two terms: a stimulated term, 
Pn (t), caused by the free evolution of Tp(x,0 + ), and 
a spontaneous term, p n D ^ (t). 

After Fourier expansion of i/j(x,0 + ), the stimulated 
Fourier components of the density are given by squar- 
ing 

OO 

iP(x,t)= £ (-l) n J 2n (0/2)e mqx - m2uJ i t . (52) 

n— — oo 

to form M 



i S) {t) = 1 /2{J2n[dsm(nw q t/2)] 

+ (-1)"J 2 „ [0cos(nuV/2)]}. 



(53) 



In particular, the average stimulated density (n = 0) is 

pf(*) = 7a[l + Jb(0)]. (54) 

Each of the stimulated terms, p n (t) , is periodic in time 
with the Talbot period, T = 2tt /uj q . 

The nth spontaneous term, p„ (t), results directly 
from the decay of the nth Fourier component of the ex- 
cited state population. Since this decay is accompanied 
by an atomic recoil, each harmonic's spontaneous term 
for n ^ acquires an additional Doppler phase which 
depends on the momentum Kk r of the emitted photon, 
nhqx ■ k r t/M. The resulting inhomogeneous dephasing 
after integrating over spontaneous emission directions 
leads to a decay of the spontaneous part of the n ^ 
Fourier components on a time scale of the order of T, 



P { nJo(t»T) 



0. 



(55) 



In the spatial domain this is equivalent to a length scale 
on the order of the Talbot distance, Lt- While we could 



use the complete results of Ref. 28 to analyze focus- 
ing including this spontaneous term for times t < T, we 
choose instead to simplify the problem by waiting exactly 
j Talbot periods (t ~ jT ^> T) , where j is a positive inte- 
ger greater than one. As a result, one can neglect joS (t) 
inEqs. @. 

The n = spontaneous Fourier component does not 
decay. For n = it follows from the Fourier expansion of 
tp e {x,0 + ), for example, that 



fh) 



(D 



} (t)= 7a [1- Ja(9)\ 



(56) 



Thus, in accordance with the conservation of probability, 
the total average density is one by summing Eqs. ( [54] ) 
and (15? 



Po(t) = pi s) (t)+pl ) D) (t) = l. (57) 
Combining these results, the total density is 

(58) 



p(M»T) = l + 5>(f> (t), 



,inqx 



We can now use Eqs. (|5l|), (|53|), (pq), and ( |58| ) as ex- 
act, Method 1 expressions to compare to an approximate 
diffraction theory. 

To present a diffraction theory of resonant focusing, 
we consider only the stimulated ground state wave func- 
tion ip (x, t) evolving from ^{x^jT + + ) = ip(x, + ) [Eq. 
Q49|)], remembering that it is first necessary to add an 
integer number of Talbot periods, jT, to reach an "ini- 
tial" time when the modulated spontaneous contribution 
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can be neglected. Moreover, we must add on the ad- 
ditional background density from the spontaneous term, 
Eq. (|56|), to the final result. 

The wave function, ip(x,0 + ) = cos[(0/2) cos (for)], is 
amplitude modulated; however, ip (x,0 + ) superposes two 
terms, each of which is phase modulated, 

ip (x, 0+) = ip+ (x, 0+) + ip- (a;, 0+) , (59a) 
ip±(x,0 + ) = V2exp[±i (6/2) cos (kx)]. (59b) 

These states can be seen as evolving independently, lead- 
ing to the components of ip {x, t) we denote as ip± (x,t) . 
Therefore, the problem of focusing after this specific type 
of amplitude modulation can be mapped onto the prob- 
lem of focusing by a phase grating. For times uo^t ~ _1 , 
> 0, near the even antinodes, kx = 0, ±2ir, . . . , the 
component i/j + (x, t) is responsible for focusing while the 
component ip- (x, t) evolves smoothly as a background 
density. Near the odd antinodes, kx = ±tt, ±37r, . . . , 
the role of the components ip±(x,t) is reversed. The 
stimulated density is then given by 



(x,t) = |V+ (x,t)+iP- (x,t)f 



(60) 



and the total density, adding on the spontaneous back- 
ground p { D) (t), by 



p(x,t) = 1 / 2 [l-J (6)]+p^(x,t). 



(61) 



Consider the atom density near x = (the even 
antinodes) for the development of the asymptotic the- 
ory. Since the period of the component wave functions 
ip± (x, t) is twice as large as the far detuned case [com- 
pare Eqs. ( |9a| ) and (59b)], the dimensionless coordinate 
and time become 



t -> LO k t, 



(62) 



noting that q — 2k and uo q = itOk- We again assume that 

1 12 3> 1. For the ip + ( x :t) component, free evolution in 
the coordinate representation including the lowest-order 
spherical aberration gives 

V>+ (x, t) = (3/6») 1/4 (4mt)~ 1,2 exp (£0/2) / (x, u) , (63) 

where the function / and variables Co and x are given by 
Eqs. (p8b|- |28dj). 

To calculate ip- (x,t), it is sufficient to consider only 
the harmonic part of the potential near x = 0, just as we 
did for the background term ([l3|) in the phase modula- 
tion case. Replacing ip(x',t') in the Fresnel-Kirchhoff 
equation (25) by ip- (x', + ) and using the expansion 
cos (x 1 ) ~ 1 — x' 2 /2, one finds 



ip. (x, t) = Va (1 + 0t) 1/2 exp (-i0/2 + ix 2 _) 
x- = [At{l + 8t)y lh x. 



(64a) 
(64b) 



This contribution remains of order unity at the focus, i.e., 
it has relative weight 9~ 1 1 i compared to Eq. (63). Since 



Eq. ( |63| ) is of higher accuracy than Eq. (64a) (relative 

corrections to Eq. ( p3|) are of order 9~ / 2 ), it is valid to 
add Eqs. (|6^) and ( |64a| ) to form the total wave function 

and calculate the density to absolute order 9 / 4 . 

Interference between ip + (x, t) and ip- (x, t) leads to a 
new effect, a Rabi-like oscillation of the focused atom dis- 
tribution. Squaring ip(x, t) — ip + (x,t)+ip- (x,t) to form 
(x ,t), we find that the asymptotic stimulated atom 
density near x = is 

p^(x,t) ~(4nt)- 1 (3/9) 1/ - 2 {F(x,Co) 

+6» 1/4 [2i/(l + 0i)] 1/a 

x [cos (9 - x 2 _) f c (x,Co) 

- sia(9-x 2 _)f s (x,Co)]}, (65a) 



f c (x, Co) + if s (x, Co) = (2n/i) 



f(x,Co) 



(65b) 



Due to the symmetry between ip + and ip—, this is the 
density near the focus at both the even and odd antin- 
odes, giving a total density of spatial period 2n/q = A/2. 
Equations ( |65| ) are the Method 2 result for the resonant 
lens. After accounting for the relevant coordinate and 
time scales (62), the dominant term, |?/>+| 2 , looks exactly 
like the far detuned result of Sec. II and leads to the same 
lens parameters if taken alone. The interference term has 
a relative amplitude 0~ / 4 near the focus (t ~ _1 ) when 
compared to the \ip+\ 2 term and oscillates sinusoidally 
with 9. 

Now we employ Method 3, the limit of this asymptotic 
result for t ~ 9~ l . Considering the interference term as 
a small correction to |j4±| 2 , we apply the far detuned re- 
sults [Eqs. © and (§3), i,a -2.20, and x f w 0.979] 
and Taylor expand Eqs. ( pq ) around Cof and Xf. The 
following asymptotic expressions can be derived for the 
focal parameters: 

t f ~ 9- 1 {i -co(wy 1/2 + 3~ 1/2 <r 3/4 

x [d 2 F{Q,Co) /dCo 2 ]' 1 [cos(0) 
xdf c (0,Q) /dCb 

- mn(e)df s (0,u,)/dw]} o=a/ 

~ _1 |l + 1.270" 1/2 +0" 3/4 
x [-0.465 cos (0) 

+ 0.110 sin (0)]}, (66a) 
p (s) (0, t f ) ~ 3 1/a [Airy 1 9^ 2 {F (0,Q f ) 

+0- 1 / 4 {[cos(0)/ c (O,^ / ) 

- sin(0)/ s (O,£/)]} 

~ O.9570 1/a |l + 0" 1/4 [0.166 cos(0) 

+ 0.703 sin (0)]}, (66b) 

w ~ 3~ 1/4 0~ 3/4 {i 
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+6- 1/4 [OF (x,u>) /dxf 1 

x [cos(0)(72/c(O,S)- f c (x,u>) 

+ (df c (0,Q)/du) {OF {x, w)/dw) 

x (d 2 F(o 1 u)/dx 2 y 1 ) 

— (the same with replacements 
cos (9)^ sin (0), fc^f.)]}*^,,^, 

0.7446*" 3/4 jl + (T 1/4 [0.198 cos (6<) 



0.162 sin (6<)]} 



(66c) 
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FIG. 4. Resonant focusing. Comparison between exact 
densities by the Fourier method (Method 1, solid lines) and 
approximate asymptotic densities near x — (Method 2, 
dashed lines) for different values of the pulse area 9. (a) Time 
evolution of the atom density at x — 0, p^ s '(Q,t). (b) Den- 
sity profile at the focal plane, p^ s \x,tf). The atom density 
is plotted excluding the constant background term caused by 
spontaneous emission, Pq (t) = (1 — Jo(0))/2. 



In Fig. ||the exact time dependence of the peak stimu- 
lated density, p^ (x, t), and the stimulated density pro- 
file at the focal plane, p^ (x, t f), a s calculated numeri- 
cally by Method 1 using Eqs. ( |51b| ) and (|S3| ) for several 
pulse areas, are compared with th e ap proximate diffrac- 
tion expression of Method 2, Eq. ( |65a|) . A comparison of 
the focal parameters as functions of the pulse area for the 
three different methods is plotted in Fig. |[ clearly show- 
ing the oscillation of the exact and approximate results 
with 6. While a diffraction theory may be unnecessary 
since an exact result can be calculated by the Fourier 
method, the physical origin of the Rabi-like term, which 



is masked by the Fourier result, has been revealed by the 
diffraction theory to be the interference of the ip+ and 
ip— components of the wave function. Further discussion 
is again reserved for Sec. V. 
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FIG. 5. Resonant focusing. Pulse area dependences of the 
focusing parameters: Peak focal densities [p(0, t/) (curves 
p)], spot sizes w (HWHM of the density profiles at the fo- 
cal planes), and focal distances tf. Computational meth- 
ods: 1. Exact Fourier (solid), 2. Asymptotic diffraction near 
x = (dashed), 3. Asymptotic diffraction near x = and 
t — (dot-dashed). The atom density is plotted excluding 
the constant background term caused by spontaneous emis- 
sion, p { °\t) = (1- J o (0))/2. 



IV. CHROMATIC ABERRATION AND 
TRANSVERSE VELOCITY DISTRIBUTIONS 

A. Chromatic aberration 

Chromatic aberration, the dispersion or wave length 
dependence of a lens' properties, results from a finite dis- 
tribution of longitudinal velocities (i.e., de Broglie wave 
lengths) in an atom beam. In trap experiments this type 
of aberration is avoided. However, for beam experiments 
detected in the laboratory frame, the atomic: distribution 
would be probed or deposited at a certain distance from 
the interaction region, L = Ut. As a result, the total 
atom density at L will be an average of each velocity's 
density at L over the flux of atoms with that velocity. 
In general, to calculate the density, we must know the 
longitudinal velocity distribution or at least its statis- 
tical properties. Previous theoretical and experimental 
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publications that have discussed chromatic aberration 
in the context of atom optical lenses have used heuris- 
tic, numerical, or Monte Carlo simulation approaches 
We attack the problem from a different 
perspective, based on averaging the exact Fourier com- 
ponents over the atomic flux distribution. The technique 
is easy to apply, the results are simple to understand, 
and the conclusions have a physical interpretation. 

For the remainder of the paper, the ideal case will re- 
fer to the monovelocity atomic beam with longitudinal 
velocity U and infinitely narrow longitudinal and trans- 
verse velocity distributions. In the previous sections we 
have presented a full characterization of the thin SW lens 
acting on the ideal beam. Two types of beams with more 
realistic longitudinal distributions are typical and will 
be considered here, (1) a velocity narrowed beam cen- 
tered around some average velocity Uq or (2) a thermal 
beam with average speed Uo- In this section we con- 
sider only focusing by the far-detuned standing wave. 
For this ty pe o f lens, the exact expression for the den- 
sity, Eq. (21t), can be averaged over the proper flux 
distribution function for a specified L by numerical inte- 
gration, provided the paraxial approximation still holds: 
Arfs,o = 2ttTi/ (MUo) <C A. In addition, the narrow veloc- 
ity distribution case allows for an approximate analytical 
solution to compare to the numerically integrated result: 
by expanding the lens parameters around the average val- 
ues of the flux distribution, we can account for chromatic 
aberration analytically using the density of the Fourier 
method. 

In the laboratory frame the atoms interact with a SW 
field with a fixed width a z along the z direction and are 
then detected at a fixed distance L from the lens. As a 
result, the pulse area 9 and the time of free flight t after 
the lens depend on the atomic longitudinal velocity U as 
U . The atomic density is defined as an average over 
the single-particle atomic flux distribution W (U). This 
distribution is a measure of the probability, W (U) dU, to 
find an atom with velocity between U and U + dU pass- 
ing through a plane perpendicular to the z axis. (One 
would divide W (U) by the atomic spatial density to get 
a true flux distribution in inverse units of particles per 
unit time per unit area.) To derive quantitative results, 
we need to specify W (U) for a given experiment. 

In this section we denote flux averages () by 



(F(U)) = / dUF (U) W (U) . 



(67) 



taking into account that only atoms having positive ve- 
locity, U > 0, have to be included. Using this definition, 
the average velocity Uo, the (dimensionless) relative ve- 
locity it, and the (dimensionless) flux width v are given, 
respectively, as 



Uo 



{U), 

(U - Uo) /Uo, 

7» 



v=[2(u 2 )] 



(68a) 
(68b) 

(68c) 



The pulse area 9 and the time of flight t for each veloc- 
ity class U can be defined by reference to the average 
velocity Uo, 

9 = 9 Uo/U = O /(1 + u) and t = t U /U = i„/(l + u), 

(69) 

where 6 ~ (-4 |x| 2 /A)(a z /U ) and t = L/U would 
be computed using Uq. The exact, flux -aver aged density 
p(x,to) follows immediately from Eq. (21b), 



p(x,t ) = (p(x,t)) 



E 



<■-' I dUW(U) 
o 



xJ n [(9 Q U /U) sin (ntoUo/U)}. 



(70) 



The integrals can be evaluated numerically, term-by- 
term, giving the flux-averaged Fourier components of the 
density. If the distribution W (U) has a non-zero width 
v, then the focusing by the SW lens is degraded. We 
show this graphically below. 

We start with the narrow velocity distribution case, 



v < 1, 



(71) 



and derive approximate, analytical expressions from Eq. 
( [70| ) to compare to an exact numerical integration. Our 
goal is to determine the effect of a flux width v on the 
focus when compared to the ideal beam results (v — > 0). 
Since v <C 1, we can expand the field area and time, Eq. 
©, as 



h (l - u + u 2 ) , t w t (l - u + u 2 ) 



(72) 



Substituting values (J72J) into p(x,t) of Eq. (21b), ex- 
panding in powers of u, and applying Eqs. (|7^) and (|6"§|), 
one finds to order v 2 that 



p [x, to) ~ p (x, t Q ) + v 2 p 1 [x, to) ■ 



(73a) 



Pi 



(x,t) = x /i ^2 {°o 4nt cos (nt) + (2 - (nt) 



n— — 00 



x sin (nt)] J' n [9 sin (nt)] + 9 2 [sin (nt) 
+ nt cos (nt)} 2 J,;' [0 O sin (nt)} } e m; 



(73b) 



where p\ (x, to) is the lowest-order correction to the den- 
sity at L as a result of the flux distribution. 

To order v 2 , four of the focal parameters (the focal 
distance tf, defined as the first maximum of the function 
p (0, to) , the peak density p(0,tf), the spot size w, defined 
as the lowest root of the equation p(w,tj) — p(0,tf)/2, 
and the depth of focus At, defined by p(0,t±) = 
p(°7/)/ 2 )> are S iven by 
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*/ — *0,/ + v *1) 



il = - 



9pi (o,t) /a* 



5 2 p (o,t) /at 2 

p(P,t f ) ~p(0,t 0J )+v 2 p 1 (0,t j); 
w ~ Wo + trtui, 

tfli = {[flp^tj/fto] -1 [VaPi (0.*) 
-Pi («>,*) 

- hd P (w,t) /dt}} t=toj w=wo , 

At = At + v 2 Ati, 

VaPi (0,to,/)-pi (0,td,+) 



[Sp(0,t)/9i] t=to + 
V2P1 (0,t ,/)-pi (0,to,_) 



[dp(0,t)/dt\ 



(74a) 
(74b) 

(74c) 
(74d) 



(74e) 
(74f) 



(74g) 



tion, 



t=tn- 



where to,/, u>o, and Ato = to,+ — to,- are the focal dis- 
tance, focal spot size, and the depth of focus, respectively, 
for #0 as calculated in Sec. II for the monovelocity beam. 
It is significant that Eqs. ( |73| ) and (frij) depend only on 
the dimensionless flux width v and not on the exact form 
of the distribution W(U). The pulse area dependences 
of the correction coefficients, t\, p\ (0, to,/), and Wi, are 
shown in Fig. |[ To verify Eqs. (73) and (|74|) for small 
v, we can now quantitatively compare the approximate 
expressions to the exact, numerically averaged density, 
Eq. ( f70| ) , for a physically reasonable and mathematically 
convenient distribution function. 




PULSE AREA 

FIG. 6. Chromatic aberration for focusing atomic beams. 
Correction coefficients, given in Eqs. J74| ) , as functions of the 
pulse area #0, assuming a narrow longitudinal flux distribu- 
tion, i) < 1: Peak focal density pi(to,/), focal distance ti, 
and spot size wi corrections. These coefficients multiply v 2 
to form the lowest-order corrections to the lens parameters. 

For example, we can find the flux distribution Wim(U) 
corresponding to the local Maxwellian velocity distribu- 



Wm> (U) = 



exp 



- [U-Uo\ /V 



(75) 



The flux distribution, average velocity, and (dimension- 
less) flux width that follow from Wim' (U) can be written 

with accuracy exp (^—U 2 /V 2 ^ 



as 



W m (U) = 



u 



■ exp 



U = Uo(l + V 2 /2U 2 ), 



U-Uo) /V 



y ^l-V 2 /2U 2 
U l + V 2 /2U 2 



(76a) 
(76b) 

(76c) 



As long as V/U g -C 1 , we can be assured that u«l. Dis- 
tribution (76a) can be inserted into Eq. ( |70| ) for p (x, to), 
and the exact density can be evaluated at the time to,/ 
for different values of the distribution width v using the 
parametrization of Eq. ( [76cj ) . Alternatively, the pre- 
cise, flux-averaged focal time, f/ = Lf/Uo 7^ to,/, can be 
found numerically from the first maximum of p(0,to) in 
time and used to evaluate the density and lens parame- 
ters. The peak reduction and spot size increase for two 
pulses areas at to,/ are plotted as a function of v in Fig. 

In Fig. ^ the density profile at the true focus, p(x,tf), 
for one of these pulse areas is shown with and without 
chromatic aberration. Both figures are explained in de- 
tail in Sec. V. 




FIG. 7. Narrow longitudinal flux distribution in atomic 
focusing. Exact reduction in the peak focal density, 
p(0, toj)/p(0, to,/) (boxes), and exact increase in the spot 
size, w/wo (triangles), as functions of the longitudinal flux 
width v for the flux distribution Wim(U) and for two pulse 
areas, 9o = 9.37 and 14.86. These pulse areas would produce 
spot sizes of 10 nm and 6.5 nm, respectively, in Na for v = 0. 
The dashed curve for p(0, to./)/p(0, to./) is t he a pproximate, 
analytical dependence for small v from Eq. (74c) or Fig. ^| 
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FIG. 8. Exact focal density profiles, p(x,tf), in the range 
\qx\ < 7r/2 for 8q = 14.86 with and/or without the effects 
of both chromatic aberration and angular divergence. The 
ideal atomic beam - v — 0, V x = - gives a spot size of 
qw — 0.139, or 6.5 rim in sodium (thin solid line). Pure 
chromatic aberration, Wim{U) longitudinal flux distribution 
- v = 0.194 and V2 ~ (dotted line). Pure angular divergence, 
P2 (Maxwellian) transverse velocity distribution - v = and 
V2 = 1 (dashed line) . Combined chromatic aberration and an- 
gular divergence - v = 0.194 and V2 = 1 (thick solid line). The 
transverse width, V2 = 1, corresponds to the single-photon re- 
coil limit. See Table 2 for other quantitative details. 

Next, we consider the other case of experimental inter- 
est, the Maxwellian thermal beam having a flux distri- 
bution 



W M (U) = 2U 3 V- i exp (-U 2 /V 2 ) 



(77) 



where V is now a thermal velocity, 
the distribution width v, still defined by Eq 
small, 



For a thermal flux, 
ic), is not 



v/64/9vr- 2 « 0.513. 



(78) 



On the other hand, since v < 1, we expect that the order 
of magnitude of the density's peak amplitude, the spot 
size (or spatial resolution) of the lens, and the focal po- 
sition should be the same. In analogy with the analysis 
above, we consider the atomic density at a distance L 
from the SW lens, where L = Vt is now defined by the 
free evolution time to for an atom with thermal veloc- 
ity V. (The average thermal flux velocity Uo is given by 
3v^FV/4 w 1.33V.) 

We can now use the thermal flux distribution (|77| ) 
in Eq. ( f70| ) to calculate the various, exact lens pa- 
rameters. To find the focus for this distribution, in- 
stead of monitoring the peak atomic density p(0, to), we 
monitor the thermally- averaged density contrast c(to) = 
p(0,t o ) /p(n,t ) from Eq (p4|). An example of c(i ) for 
do = 23.0 is shown in Fig. [|. The focal plane posi- 
tion is defined by the time when c(to) is largest, to = tf 
(L = Vtf). The atomic density at tf, p(x,tf), for this 
pulse area is also shown in Fig. |9|. Several lens param- 
eters as functions of 9q are shown in Fig. |ol including 



the focal time tf, the contrast at the focal plane c(tf), 
and the spot size w. 
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FIG. 9. Thermal longitudinal 

beam focusing, (a) Time evolution of the atomic density con- 
trast, c(to) = p(0, to)/p(ji, to), and the density profile at the 
focal plane, p(0, tf). (b) Atomic density profile at the focus, 
p(x,tf), plotted in the range — A/4 < x < A/4. The pulse area 
for these curves, 9o = 23.0, is needed to produce a spot size 
of w — 6.5 nm in Na using a thermal beam. 




PULSE AREA 

FIG. 10. Thermal longitudinal beam focusing. Exact pulse 
area dependence of the focusing parameters: Focal density 
contrast c(t/) = p(0, i/)/p(7r, tf), focal distance tf, and spot 
size w (HWHM of the density profile at the focal plane). 
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B. Transverse velocity distributions: Angular beam 
divergence/ Trap velocity distributions 



Transverse velocity distributions put an additional 
constraint on the thin SW lens performance and have 
been a limiting factor in Fresnel atom optics and inter- 
ferometry. For an atom beam, the angular divergence 
is equivalent to the transverse velocity distribution of 
atoms in a cold trap, again allowing a time domain treat- 
ment. The transverse and longitudinal beam velocities 
are generally related linearly by the angle in the parax- 
ial limit. This is the case both for atomic beams colli- 
mated by apertures and for beams that are laser-cooled 
in the transverse direction. Given that we have already 
discussed chromatic aberration in detail, here we present 
comparative focal parameters and densities which include 
either transverse effects alone or simultaneous longitudi- 
nal and transverse velocity effects. Only the narrow lon- 
gitudinal velocity atom beam will be considered in detail 
since this type of beam leads to a better focus. 

We will look at two cases. For both cases the Fourier 
method suggests a straightforward interpretation of the 
effects of a transverse velocity distribution in terms of an 
inhomogeneous decay of the density spatial harmonics. 
First, we treat the case where the beam is monochro- 
matic as in Sec. II and the transverse velocity distribu- 
tion is either a uniform or one-dimensional Maxwellian 
(thermal) distribution of transverse velocities. The uni- 
form distribution might apply to an atomic beam colli- 
mated by a slit or a pair of slits while the thermal distri- 
bution might apply to a laser-cooled beam or to atoms 
released from a laser-cooled trap. The resulting focal pa- 
rameters and densities for the slit-collimated and laser- 
cooled atoms are similar qualitatively. Then, we consider 
a more general beam case where the transverse velocity 
and longitudinal flux distributions are decoupled, but the 
longitudinal velocity dependence of the time couples the 
transverse integration to the longitudinal. This applies 
to atomic beams that are laser cooled before interacting 
with the SW lens. The theoretical treatment of each of 
these cases is facilitated by the Fourier method. 

Our first case is the monovelocity longitudinal beam 
(v — > 0) with a transverse velocity distribution P(U X ). 
To unify the longitudinal and transverse results, we start 
with Eq. @ and take W(U) = S(U - U ). This is 
equivalent to starting from Eq. (21a) or ( |21b| ) and tak- 
ing U = Uq, t = to, and 6 — 8q. Assuming we are 
describing the atomic beam rather than the trap, the 



oo 

>(x,t )= ]T e ma; J„[#osinK,)] J dU x P(U x ) 



atoms are again detected, deposited, or used for lithog- 
raphy at L — Uat(). The inclusion of an initial transverse 
velocity U x during the atom-field interaction will lead to 
a Doppler shift of the field frequency. (Note that the 
dimensionless transverse velocity scale that is consistent 
with the dimensionless variables (20) is the single-photon 
recoil velocity, Vk — hq/2M.) We rigorously account for 
the Doppler shift by taking x — ► x + U x t and averaging 
the density over P(U X ), 



(79) 



The assumption is that the initial transverse density ma- 
trix is diagonal in momentum space and can be described 
by a normalized velocity distribution P(U X ). In terms of 
the longitudinal velocity Uq of the beam, the transverse 
velocity can be written as U x — Uqip, where if is the 
propagation angle with respect to the z-axis. 

The two transverse distributions of interest, the uni- 
form distribution P\{U X ) and Maxwellian distribution 
P2{U X ), take the following form for V x <C Uq: 

1. P 1 (U x ) = {2V X )~ 1 for U x € [-V x , V x ], (80a) 

2. P 2 (U X ) = (v^Kr 1 eM~(U x /V x ) 2 ]. (80b) 

Each can be written as a distribution over beam angles <p, 
if desired, by defining the divergence angle <fd = V x /Uq. 
The first distribution might be formed when a thermal 
transverse distribution, peaked at U x = (ip = 0) with 
velocity width much greater than Vk, is collimated by a 
pair of slits which select a small range of propagation an- 
gles up to ifd- This beam has a rms transverse velocity 
spread of V\ — V x /y/3. The second distribution might 
be formed by transverse laser cooling, with or without 
previous slit collimation, and has a rms velocity width 
of V<i — V x /y2. (Technically, the laser cooling process 
can create a non-diagonal density matrix with transverse 
momentum state coherences that are not accounted for 
here but can be incorporated into a more general Fourier 
result. We assume the distribution is narrow and diago- 
nal.) 

Fortuitously, both distributions again give closed-form 
Fourier coefficients for the density. Inserting Eqs. ( |80| ) 
into Eq. (|79|), one finds 

c sm[nV x t ] 
nV x t 

x J n [0 O sin (raf )] ; (81a) 

OO 

2. p(x,to) = exp[ma; - (nV x t /2) 2 ] 

n— — oo 

xJ n [9 Q sin (nt )}. (81b) 

The interpretation is the same as for an atomic free in- 
duction decay experiment fUjfl . The harmonics of the 
total density undergo an inhomogeneous decay as each 
velocity class U x evolves with its own Doppler-shifted 
frequency. The Pi result was recently explained in de- 
tail as it is isomorphic to the evolution of a model, one- 
dimensional BEC or a degenerate Fermi gas after inter- 
acting with a SW pulse |^{| . A similar expression to the 
Pi result was used recently to explain the decay of peri- 
odic echoes |l3| and the decay of a quasiperiodic atomic 
focusing and Talbot scheme p7fl . For both distributions 



i. p(x,t Q ) = ^2 
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the harmonic decay functions of Eqs. ( |81| ) depend on 
nVxOq 1 near the focus, to ~ ^q' 1 . For V x % 1 only har- 
monics greater than n max ~ 8 undergo significant damp- 
ing, so it appears that the transverse velocity width can 
be up to the order of the recoil velocity to avoid a signif- 
icant decay of the high-order Fourier components neces- 
sary for a highly-peaked focus. However, the restriction 
V x % 1 may be sufficient but not necessary. We have 
found in the previous sections that the lens parameters 
scale with 6 , and, therefore, in Sec. V we show that the 
restriction on V T is somewhat relaxed. 



Th e foc al di stan ce t f , defined as the first maximum of 



Eq 



la|) or ( plb| ), moves closer to the lens for V x > 
when compared to the focal distance for V x = 0. This 
makes sense since the transverse distribution causes a 
decay in the modulated density in time. However, this 
shift of the focal distance is secondary in importance to 
the changes in the peak density at the focus and the spot 
size. Moreover, as 8 increases, the focal distance includ- 
ing the transverse velocity spread converges rapidly to 
the ideal value near to ~ O^ 1 . We have calculated the 
peak density and spot size for selected values of V x and 
compared them to the ideal case. In Fig. [ll] the peak 
focal density, p (0, tf) , and spot size w are plotted versus 
#o for V x = (ideal), V\ = Va = 1 (recoil limit), and 
Vi = V<z = 2 (twice the recoil limit). For each of the 
five curves, the exact, different focal time at each value 
of 8q is used to calculate the parameters. For the same 
rms velocity spread, the Pi distribution gives a narrower, 
more peaked density spot than the P\ distribution. It is 
significant that the percentage error in the peak density 
and spot size is decreasing as 9q increases. This result is 
expected since inhomogeneous decay is less pronounced 
at earlier times, to ~ 67, 1 . 



To treat the full problem of focusing an atomic beam 
with velocity distributions more generally, we must ac- 
count for the angular relation between the longitudinal 
and transverse velocities. The atomic velocity in a beam, 
U = Uz + U x x, can be written U ~ |U| z + |U| cpx 
for small angles ip, where |U| ~ U to lowest order in 
(p in the paraxial approximation. Since atomic beams 
are inherently formed by a process which truncates ei- 
ther the divergence angle or the transverse velocity dis- 
tribution, it is not a simple matter to describe them 
theoretically. However, if the beam is localized around 
some large longitudinal speed and then laser cooled in 
the transverse direction before the SW lens, the longi- 
tudinal and transverse distributions are effectively in- 
dependent. We can model the total (normalized) flux 
probability W(U)dU as a product of the two indepen- 
dent distributions, WiM(U)P2(U x )dUdU x , using the lo- 
cal Maxwellian longitudinal flux distribution Wim(U), 
Eq. (76a), a nd th ermal transverse velocity distribution 
P 2 (C4]7Eq. 
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PULSE AREA 

FIG. 11. Angular divergence (transverse velocity distribu- 
tion) in atom focusing, neglecting chromatic aberration. Ex- 
act pulse area dependence of the lens parameters: Peak focal 
density p(0,tj) and spot size w (HWHM of the density profile 
at the focal plane). The solid curves are for the ideal atom 
beam (with no angular divergence), the dotted curves are 
for the uniform transverse distribution (Pi), and the dashed 
curves are for the Maxwellian distribution (Pn). The middle 
pairs of curves, closest to the ideal case, are for Vi = V% = 1, 
the rms single-photon recoil limit. The outer pairs are for 
Vx = V2 = 2, twice the recoil limit. 

Since the time of flight, t = L/U, to the detection 
plane, L = Uoto, depends on the longitudinal velocity, 
the Doppler phase e mUxt ° Uo ' u does as well. Inserting 
the integration of the Doppler phase over the distribu- 
tion PziUx) into Eq. (|7p|), we find the velocity- averaged 
density to be 

00 />oo 

p{x,t )= Y, dUW(U)exp[~(nV x t Uo/2U) 2 } 

x { J n [{0 o U a /U) sin {nt U a /U)}} . (82a) 

The Fourier component of the density now has an addi- 
tional Gaussian multiplier to integrate over longitudinal 
velocities. If the longitudinal width v is small but non- 
negligible and we want V x < 1 (laser cooling to the recoil 
limit), the argument of the Gaussian is on the order of 
the argument of the sine function. As a result, we cannot 
neglect the Gaussian decay in the integral by assuming 
W (U) is sharply peaked around Uq and thus evaluating 
it at U — Uq (v — > limit). In Sec. V we treat two 
examples of experimental interest with this expression, 
an atomic beam which is laser-cooled to the transverse 
recoil limit, both with and without a longitudinal flux 
width v. This is shown in Fig. || along with the ideal and 
purely longitudinally-broadened cases. 
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V. DISCUSSION 

In this article we have developed a theory of atom 
focusing by standing wave light fields in the thin lens 
regime. We have shown that the exact analytical expres- 
sions for the Fourier components of the atomic density, 
recently used in the theory of atom interferometry, are 
useful for numerical calculations of the focusing effect 
and all its relevant lens parameters, including the effects 
of chromatic aberration and angular divergence. We call 
the Fourier technique Method 1. 

Thin lens focusing becomes especially effective for a 
large field area 6. When y/0 ^> 1 for the monoveloc- 
ity case, using the Fresnel-Kirchhoff diffraction integral, 
we can consider only the first anharmonic term of the 
phase shift created by the SW light field when expanding 
around the potential wells at x' — lixm. This procedure, 
labeled Method 2, leads to a manageable expression j30| ) 
for the atomic density distribution which can be used to 
evaluate the focal parameters by numerical integration. 

Within a small vicinity of the focal time, Eq. (|30| ) 
simplifies further to the asymptotic density profile (|36[). 
This is Method 3. The analytical correction to the focal 
time and analytical expressions for the peak focal den- 
sity, spot size, and depth of focus as functions of the 
pulse area 8 follow from this result. The last three lens 
parameters scale as (9 1 / 2 , 6> -3 / 4 , and # -3 / 2 respectively, 

while the focal time scales as B^ 1 ^1 + 1.270 - Z 2 ^ [see 

Eqs. @, @, (|l]), and @]. From Fig. | one sees 
that in the range < 8 < 40, the accuracy of the Method 
3 asymptotic expressions for the peak focal density, focal 
time, and spot size can be as small as 17%, 1.9% and 
18% respectively. The asymptotic parameters obtained 
by Method 2 using Eqs. (^0|) are generally more pre- 
cise, giving accuracies for the same parameters as small 
as 5.8%, 0.24% and 4%, respectively. 

The accuracy of the asymptotic result for the focal 
time tf is an order of magnitude better than for other 
focusing parameters. Evidently, this is a consequence of 
the fact that a finite value for the focal time, t ~ 0~ l , 
comes out of the harmonic approximation for the lens 
potential. The first anharmonic term then leads to a 
relative correction of the order of O^ 1 / 2 [see Eq. (5^)], 
and higher-order anharmonic corrections (not included 
in our consideration) should have relative weights # _1 
and larger inverse powers of 6. As for the other param- 
eters, which can be derived only by including the first 
anharmonic term, the next order corrections should have 
a larger relative weight of Q~ x l 2 . 

For example, the constant offsets of the asymptotic 
peak densities from the exact result in Fig. ^ can be ex- 
plained as follows. A correction to the asymptotic wave 
function at x = [Eq. (31)] of relative weight O^ 1 / 2 
would result from including the next order spherical aber- 
ration, the x' 6 term in the expansion of cos(x') in Eq. 
(28a). When the wave function is squared to form the 
density, this relative correction of weight Q^ 1 / 2 multiplies 



the lowest-order peak density that grows as # 1//2 to give 
a constant term, which is automatically accounted for 
by the exact Fourier method. Hence, the constant off- 
set of the asymptotic peak density is not included in our 
asymptotic expressions. 

In this article we have also analyzed a new type of 
atom focusing which arises when the SW field is resonant 
with the atomic transition. The density profiles and fo- 
cal parameters for the resonant lens are shown in Figs. 
|] and ||, respectively. To gain insight into the different 
interactions that can lead to focusing, we can consider 
the far-detuned and resonant cases from a more general 
point of view than that taken in Sees. II and III, still as- 
suming that the field envelope is short enough to neglect 
spontaneous emission during the pulse. When an atom 
interacts with a pulsed light field, the states of the sys- 
tem can be decomposed into a set of semiclassical dressed 
states of the atom plus field. If the atom-field interaction 
is adiabatic or the pulse turns on instantaneously, these 
dressed states are instantaneous eigenstates of the total 
Hamiltonian and therefore undergo a phase modulation. 
When the light field modes have a modulated intensity, 
the dressed state energies and, therefore, phase evolu- 
tions are also spatially modulated. Focusing can occur 
near the intensity extrema, which correspond to spatial 
minima of the dressed state energies. If initially the bare 
atom was in a pure state such as the ground state \g), the 
atomic wave function after the interaction is generally a 
coherent admixture of the two dressed states, which have 
an energy separation that is spatially modulated. There- 
fore, an interference term in squaring the bare state am- 
plitudes can lead to a Rabi-like oscillation of the total 
density and its pr oper ties. This Rabi-like oscillation can 
be seen in Eq. (35a) and Fig. [| for a resonant SW 
field. Nonadiabatic, detuned atom-field couplings can 
lead to a similar effect. The evolution of the pure and 
dressed states in this regime has been considered previ- 
ously without center-of-mass spatial effects 1 43 . For the 
far detuned standing wave (|A| ^> \x\) or an adiabatic 
turn-on of the field, the ground state evolves into only 
one of the dressed states while the other dressed state 
has a negligible amplitude. As a result, one does not 
observe any interference or field dependent oscillation in 
this case. 

Up to this point in the discussion, we have reviewed the 
ideal situation of the ideal atomic beam having no angu- 
lar divergence and no longitudinal velocity distribution. 
Our calculations in the thin lens regime have shown that 
focusing results in a relatively large peak focal density. 
For the most part, this result is in contrast to previous 
experiments carried out in the thick lens regime, where it 
proved difficult to achieve such large density peaks. The 
exceptions are certain experiments in Na at A = 589 nm 
on the 351/2, F = 2 -> 3P 3/2 , F' = 3 transition In 
these thick lens experiments for a thermal beam having 
thermal velocity V ~ 8.6 x 10 4 cm/s, the spot size reso- 
lution and contrast were w = 10 nm, c(t f ) = 10 @ and 
w = 6.5nm, c(tf) = 6 0, respectively. As an exercise, 
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we can compare that experimental data using a thermal 
beam and thick SW lens with our results using a ideal 
atomic beam and thin lens. From Fig. || or an exact 
calculation, a pulse area of 9 ~ 14.86 (9.37) is needed for 
a resolution of qw = 0.139 (qw — 0.139), or w — 6.5 nm 
(10 nm) in sodium. If we set the beam velocity U equal 
to 8.6 x 10 4 cm/s, then focusing for 9 ~ 14.86 occurs at a 
distance Lf = Utj = 127 um. The atomic density profile 
at this distance is shown with the thin solid line in Fig. 
|§|. The contrast at the focus, defined by Eq. (Q) as the 
ratio of focal atomic densities at qx = and qx = tt (not 
shown in Fig. |), is c(t f ) = 27.2. 

We can determine the standing wave power needed to 
achieve a pulse area of 9 = 14.86 for the two-level system. 
This system can be realized in Na using optical pumping 
to the ground state sublevel having magnetic quantum 
number m F = ±2. A circularly polarized standing wave 
drives only the 3S 1/2 , F = 2, m F = ±2 -> 3P 3/2 , F' = 
3, nip* — ±3 transition which has the dipole moment 
matrix element \p\ — 6.37 x 10~ 18 esu ■ cm. To achieve 
the intensity necessary for atomic focusing, the SW field 
must be formed by focused laser beams. In the exper- 
iments of Refs. the fields have been focused to a 
circular spot having radius a z ~ 29 pm p4| . For experi- 
ments in the thin lens regime, it is sufficient to focus the 
field only in the z-direction, along which the atomic beam 
propagates. We assume that the laser field intensity has 
a Gaussian profile in the (y, z) plane, 



\E (y, z)\ 2 = \E\ 2 exp (-2y 2 /a 2 y - 2z 2 /a 2 ) 



(83) 



where a x and a y are field radii along the z- and y-axes, 
respectively (i.e., the directions perpendicular to the laser 
beam propagation along the x axis). For homogeneous 
atomic focusing into a set of lines, one should choose a y 
to be larger than the atomic beam radius. From the 
data in Refs. o\, should be as large as 0.5 cm. The 

general field mode ( p3| ) can be created using cylindrical 
optical lenses. The far-detuned pulse area 9 (10) can 
be re-expressed through the constituent traveling wave 



field powers of Eqs. 
as 9 = -0F2 7 / 2 |/x 



(P and (|3j), P = c\E\ z a z a y /16, 
2 P/(h 2 ca y AU). For further esti- 



mates we assume that a z = 29 pm, a y — 0.5 cm, and 
A = -27T- (1.71 GHz). Then, for a field area of 9 = 14.86, 
the required laser power is P = 1.89 mW. Even this power 
level is four times less than the power used in the experi- 
ment of Ref. (7| . To compare directly to the circular focus 
of the field in Ref. M, we must take a y — 29 um, giving 
us a required power of P = 11.0 pW. These relatively low 
power levels are sufficient for our consideration of thin 
lens focusing. 

We can now discuss the influence of chromatic aber- 
ration assuming no angular divergence. A longitudinal 
velocity distribution in the atomic beam leads to a degra- 
dation of the focus. But even for the thermal beam, the 
normalized distribution width v is still less than 1 [see 
Eq. ([78])], and we have shown that the order of mag- 
nitude of the focusing parameters is the same as for the 



ideal beam (v — 0). For narrow flux distributions, v <C 1, 
the lowest-order corrections to the atom density behave 
as v 2 [see Eqs. j73])] if 9 > 1. This is a simple conse- 
quence of the fact that terms linear in the deviation from 
the average atomic flux vanish after averaging over veloc- 
ities for any symmetric, or asymmetric, distribution. In 
particular, this principle applies to the analytical Fourier 
components of the density multiplied by the narrow, lo- 
cal Maxwellian flux distribution Wim{U), Eq. (fro|), for 
small v. 

In Fig. we show the exact (using Wim(U)) and ap- 
proximate reductions of the density at the ideal beam 
focal time, p(0,toj)/p(0,to t f) [Q, as functions of v for 
two values of 9a, 14.86 and 9.37. Again, these choices for 
6 are the pulse areas required to achieve wq = 6.5 and 
10 nm spot sizes H0], respectively, using ideal sodium 
beam s. The quadratic dependence, consistent with Eqs. 
(74c), is evident for small v, but the exact density peak 



actually reduces more slowly (linearly) at larger values 
of v where the small v approximation breaks down. At 
this point the expansions of the Bessel functions of Eq. 
( fFoj ) near toj ~ 9q 1 in terms of u are invalid since 
n max > 1, or equivalently v ~ Oq 1 '- a term propor- 
tional to I it I would appear in an expansion around the 
average velocity. We take this as evidence that near this 
value of v for a fixed 9q , the focal properties go from be- 
ing dominated by spherical aberration to a regime where 
a more complicated combination of spherical and chro- 
matic aberrations is important. 

To test this idea, we have also examined the depen- 
dence of the spot size increase, w/wo, on v for 9o = 14.86 
and 9.37 using the exact density for Wim(U). In Fig. 
|t] the spot size's quadratic dependence for small v turns 
into a linear dependence near the same value of v that 
the peak density reduction deviates from the approxi- 
mate theory. We see that the chromatic aberration has 
a smaller effect on the spot size than on the peak den- 
sity. Note that from a ray tracing argument for far-field 
focusing through an effective parabolic lens aperture of 
full width ~ A/2, a strong, ^-independent, linear depen- 
dence of w ex ttv (or equivalently (\u\)) would be expected 
for a focus dominated by chromatic aberration in the ab- 
sence of spherical aberration |19|. However, this is not 
the case for the SW field lens, where spherical aberration 
determines the limiting resolution up to large values of 
9q. Thus, the chromatic aberration we find is less severe 
and depends strongly on #o- We have not explored this 
#0 dependence further. However, from Fig. ^ we deduce 
that if we fix v (and by implication the atomic beam 
properties), a threshold pulse area, 9 t n ~ must exist 
for each v to mark the breakdown of the small v expan- 
sion. For 9q ^ 9th the effects of chromatic aberration 
should be calculated exactly. 

To further demonstrate the small v results for a nar- 
row longitudinal distribution, we can define a to be the 
percentage reduction of the peak focal density as a func- 
tion of #0 and v. The flux distribution width for a given 
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9q and a is therefore v 
Eq. 



y/ap(0,t o j) /p! (0,t O j) from 
(|74c[) . From Fig. |6j or for 6 = 14.86, to limit 
the peak density reduction to a = 10%, we can esti- 
mate that a normalized flux width of v ~ 0.086 would be 
needed. (For the local Maxwellian distribution, Wim(U), 
a flux width of v ~ 0.091 is found to be sufficient by an 
exact calculation with Eq. j70|).) Using v ~ 0.086 to 
calculate the other parameters approximately from Eqs. 
( [74b| ) and (74c), the percentage corrections to the focal 
distance and spot size are even less than a, 3.9% and 
1.3%, respectively. A longitudinal width, v ~ 0.086, is 
typical of the supersonic beams produced by seeding an 
inert (noble) gas supersonic expansion with sodium fl45"[ . 
While the mean longitudinal velocity of a beam produced 
recently from a BEC of sodium atoms by Bragg scatter- 
ing 1 14 1 was only 2hk/M, or 6 cm/s, the rms longitudinal 
and transverse velocity widths were only approximately 
0.16?ifc/M and 0.30?ifc/M, respectively, giving v ~ 0.11 
and justifying the paraxial approximation. In these ex- 
periments, mean longitudinal velocities up to 11.9 Tik/M, 
or 35 cm/s, with similar widths were also achieved using 
higher-order Bragg scattering, implying v <C 0.1. Since 
experiments producing atomic beams of this type are in 
their infancy, we expect that achieving relative velocity 
widths of v <C 0.1 will shortly become routine. 

In Fig. ^ the atomic distribution at the focus p (x, tf) 
(dotted line), as calculated exactly using Wim(U) in E< 
( f70| ) for 6*o = 14.86 and our extreme case from Fig 
of v ~ 0.194 (V/Uq = 0.2), is compared to the ideal 
case (thin solid line). While the reduction of the peak 
density from 11.4 to 7.81 is a significant 31% on a scale 
where ~ 0.5 is the background at the focus, the half- 
width defining the spot size is broadened by only 5% (or 
0.33 nm in sodium), and the coherent oscillations of the 
density along x are still prevalent. This suggests that 
even for v w 0.2 and a focal peak reduced by chromatic 
aberration, a high resolution, large contrast focus is pos- 
sible with the thin SW lens. In addition, using the flux 
distribution Wim(u) with v = 0.194, we have also calcu- 
lated that the field area, 0o = 15.7, is required to restore a 
spot size of qw = 0.139, or 6.5 nm in sodium. This corre- 
sponds to the laser power, P — 2.0 mW, for a z — 29 pm, 
o y = 0.5 cm, and A = -2tt • (1.71 GHz). 

We have also considered the focusing of a thermal 
atomic beam with flux Wm(U). The results of the calcu- 
lations are graphed in Figs, p] and [h]. For this part of the 
paper only, we changed the definition of the focal distance 
to be one where the contrast (^J) , and not the density, is 
optimized. The contrast, c(to) = p (0, to) jp (ir, to) , for 
6» = 23.0 at a distance L from the SW field (t = L/V) 
is shown in Fig. ||. We see that c (to) contains three 
local maxima of approximately equal weight near the fo- 
cus. They arise as a result of the time-oscillations of the 
flux-averaged, background atom density near to = tf, 
p(ir,to). These contrast maxima "compete" with one 
another in some sense as different values of the pulse 
area correspond to different velocity classes and there- 



fore different focal distances. This effect results in the 
discontinuities in the pulse area dependences of the lens 
parameters seen in Fig. [Tc| . 

For the longitudinal thermal beam and 9o = 23.0, the 
focal contrast is equal to c(tf) = 8.03 while the spot 
size or spatial resolution is again qw — 0.139. For the 
same experimental parameters as above, a z = 29 pm, 
(T y = 0.5 cm, and A = — 2ir ■ (1.71 GHz), a laser field 
power of P = 2.9 mW is needed to produce a pulse area 
of O = 23.0. 

Finally, the atomic beam angular divergence has a cru- 
cial impact on the focus, independent of the chromatic 
aberration. For a beam with no chromatic aberration, 
U = Uo, if atoms move at an angle ip with respect 
to the z-axis, the density profile displaces a distance 
Sx ~ Uo<ptf ~ U x tf along the x axis. This displace- 
ment should ideally be much smaller than the funda- 
mental spot size w that accounts for spherical aberra- 
tion. We can use the previously generated, exact spot 
sizes and focal times for each 6*o to estimate the max- 
imum allowed rms transverse velocity from this argu- 



ment, V r 



(/) 



'j/tf. This corresponds to an angular 



divergence iff = Vrms/Uo ~ w/(tfU Q ). We expect that 
allowing V rms to be Vrms will at most double the spot 
size. 

Alternatively, from Eqs. ( |39| ) and ( |4l"| ) for dimension- 
less tf and to, the asymptotic restriction on the angular 
divergence, ip < iff, is 



iff 



0.744- 



hq 



2MU Q 



?o /4 (l - 1.270Q V2 ) 



0.744—. 

Uo 



V k 1/4 



(84) 



It is interesting to note that the asymptotic focusing re- 
striction is ~ 1 / 4 times less severe than the recoil limit, 



2MUq ~ A : 



(85) 



which is the condition required to observe atom interfer- 
ence in the single interaction re gion g eometry (for exam- 
ple, the atomic Talbot effect p^ , p5| , p3| , p4] |). Condition 
( p5J ) arises from the requirement that atoms moving at 
an angle (p not be displaced more than A/2 at the Talbot 
distance Ly ~ A 2 /2A ( jb. Even though effective focus- 
ing requires the displacement to be 3 / 4 times smaller 
asymptotically (since qw cx -3 / 4 ), this displacement oc- 
curs at a distance that is 0-times shorter than the Talbot 
distance (since Lf — Uotf ~ Lt/0o)- 

This argument is consistent with Eqs. ( |si| ) and Fig. 
|l"l| . This figure suggests that small spot sizes with sharp 
density peaks can be achieved for laser powers in the 
microwatts to milliwatts range even for transverse distri- 
bution widths larger than the recoil limit. (In Fig. [Tl| , 
note that qw = 0.3 corresponds to a 14 nm spot size in 
sodium.) For example, taking 0o = 40 and remember- 
ing that the slit-collimated beam has a rms transverse 
velocity V\ (V X V3) in units of Vk, the exact spot size 
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for Vi = is qw = 0.0577 from Eq. (glaj). Now, as- 
suming the divergence is the maximum allowed, we set 
Vx = w/t f = 1.883 for 6 = 40 and find that the spot 
size is broadened to 0.100 (4.69 nm in sodium), or 73%. 
The spot size has less than doubled for a transverse width 
almost twice the recoil limit. Similarly, the asymptotic 
predictions for 9q = 40 give an increase from a spot size 
of qw = 0.0468 (not 0.0577) for Vi = to a spot size 
of 0.0826 for Vi = U <ff = 1.495 from Eq. (f§) in Eq. 
(31a), a 76% increase. The broadening is even less severe 
for Maxwellian transverse distributions. 

We now take a case of experimental interest to demon- 
strate the combined effects of transverse and longitudi- 
nal broadening on the focal density. The pulse area, 
9q = 14.86, is used again. We show results for the Pi 
(Maxwellian) transverse distribution, Eq. ( |30b| ), since 
laser cooling has been shown to improve the focus in thick 
lens experiments In Fig. (p|) we plot the density 

at the focus p(x,tf) (dashed line) from Eq. ( jSlbj) for a 
monochromatic beam (v — > 0) cooled to the rms single- 
photon recoil limit, V% = 1 (V x — \/2). In addition, wc 
show the atomic focal density (lower solid line) from Eq. 
flB2al) for 6> = 14.86, v = 0.194 (V/U Q = 0.2 in W m {U)), 
and V2 = 1. The latter curve combines the attributes of 
the other curves in this figure to provide a real picture 
of the focus. The following table summarizes the results 
for Fig. |: 

Table 2. Focusing of an atomic beam for 6q = 14.86. 
A velocity-averaged spot of qw — 0.190 corresponds to 
8.91 nm in sodium, for which Vk — 2.95 cm/s and 
Lu q ~ 6.29 x 10 6 rad/s. 



rms velocities smaller than the recoil velocity |14j of- 
fer a promising avenue to observe and characterize thin 
lens focusing effects. We would even argue that focal 
patterns nearly identical to those described in this pa- 
per have already been achieved but not recognized in 
recent condensate experiments on the so-called momen- 
tum space Talbot effect [[l6| . In these experiments a cloud 
of cold atoms from a condensate interacts with two off- 
resonant SW pulses, separated by a variable time delay. 
The detection scheme is insensitive to Fresnel effects as 
the researchers image the far-field diffraction pattern and 
therefore the atomic momentum distribution. However, 
the accurate fit of the data in that experiment to our 
Raman-Nath theory suggests that the thin lens focusing 
effect is occurring after the first pulse along the SW grat- 
ing direction. A detection scheme sensitive to the Fresnel 
density pattern may be necessary to image the focusing 
atoms in the cloud. 
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Our work suggests emphatically that the thin SW lens 
can focus atoms effectively if the angular divergence is 
cooled near the recoil limit, even in the presence of 
strong chromatic aberration. This makes physical sense if 
one considers that the longitudinal velocity average over 
slightly different focal regions for each velocity subclass 
is a slowly-varying integral over the Fourier amplitudes 
when compared to the more sensitive transverse veloc- 
ity integral, an average over the Doppler phases. Atomic 
beams of the type considered throughout this article can 
be made in the laboratory with current technology. 

In addition to the experimental possibilities for stan- 
dard atomic beams that limit angular divergence by slit 
collimation or transverse laser cooling, recent experi- 
ments that use higher-order Bragg diffraction of a BEC 
to form atomic beams with transverse and longitudinal 
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